Discrete symmetry's chains and links between 
integrable equations 



A.V. Yurov 

Theoretical Physics Department, 
^ ! Kahningrad State University, 

O ! 236041, ALNevsky St., 14, Kahningrad, Russia 



(N 
(N 



o 



email ym'ov@freemail.ru 



Abstract 



The discrete symmetry's dressing chains of the nonlinear Schrodinger equa- 
^ ' tion (NLS) and Davey-Stewartson equations (DS) are consider. The modified NLS 

(mNLS) equation and the modified DS (niDS) equations are obtained. The explic- 
itly reversible Backlund auto-transformations for the mNLS and mDS equations are 
constructed. We demonstrate discrete symmetry's conjugate chains of the KP and 
DS models. The two-dimensional generalization of the P4 equation are obtained. 



o _ 

^ ■ 1 Introduction 

o 

I Darboux transformations (DT) are very useful method to construct exact solutions of 

integrable PDE [1]. Some times ago a new approach to the DT (called dressing chains 
of discrete symmetries) emerged at the horizon of soliton's mathematics. This approach 
! to proliferate integrable equations was proposed in the example of the Korteweg-de Vries 

^ I (KdV) and sine-Gordon equations [2]. The scheme to proliferate the integrable equation 

(we choose the KdV for definitenes) starts out from LA pair and DT of this equation. 



■ After very simple transformations (see Sec. 2) we obtain new integrable equations with 

^ their LA pairs. The modified KdV (mKdV=m^KdV), m^KdV and m'^KdV equations 

were thus constructed, the second of which can be reduced to the exponential Calogero- 
Degasperis (CD) equation by an exponential point change of variables and the third 
contains the elliptic CD equation. The only price we have to pay for this simplicity is a 
rapidly growing amount of calculation. 

Despite its apparent simplicity, the method of dressing chains is an extremely powerful 
method, as can be illustrated in following example (see details in [2-4]): as noted in [2] 
and [3], the m°KdV equations (n=0,...,3), together with the Krichever-Novikov equation, 
exhaust all the integrable equations of the form ut + u^xx + f{uxx,Ux,u) = 0. In [4] we 
generalized the dressing- chain method to considerably (l-|-2)-dimensional nonlinear equa- 
tions: the Kadomtsev-Petviashvili (KP) and the Boiti-Leon-Pempinelli (BLP) equations. 
A new result, which is characteristic precisely of multidimensional equations, is the dis- 
covery of two types of dressing chains for the KP equation (we call these the conjugate 
chains). The new chain can be constructed with the help of binary DT. 



The aim of this work is the generahzation of the dressing chains method for the 
nonUnear Schrodinger (NLS) and Davey-Stewartson (DS) equations. The plan of this 
paper is as follows. In Sections 2 and 3 we review the dresing chains method for the 
higher and lower KdV equations and for the higher KP equations. The dressing chains 
for the NLS and for the DS equations is discussed in Sections 4, 5. In Sec. 6 we develop 
a technique for constructing of exact localized solution of DS equations with a reduction 
constraint imposed and give some of these solutions. We conclude with a discussion on 
possible role of discrete symmetry's dressing chains in the theory of integrable PDE. 

And last but not least. The dressing chains method allow us to construct new inte- 
grable equations and, in the other hand, to find a links between known ones. For example, 
the sine-Gordon, KdV, and CD equations can be constructed starting out from the KdV 
equation via dressing chains. Maybe it is only one integrable equation really? 

2 2D-P4 equation 

Let us begin our discussion of the Borisov-Zykov method by analysing the simplest of all 
possible integrable systems, the KdV equation. The starting point is the LA pair for the 
KdV equation 

tpxx = (m - A)^, tpt = 2{u + 2\)ip^ - u^il^, (1) 
where the KdV equation 

Ut - 6uu^ + U^^a: = 0. (2) 

Setting r = 4>x/4>j where is a partial solution of (2) with the X — /i we can write the 
LA pair as 

rro = -r'^ + u-n, Tt = [2{u + 2h)t - Ux]^ . (3) 
Exepting u from (4) we get 

Tt = 6tVj; - T^xx + GlITa;. (4) 

li jj, = then we have the well-known mKdV equation. To construct LA pair for the (4) 
we use the DT, 

Setting (7 = ■0i,a;/'0i) we get the x-chain 

{a + r)^ = -a^ + T^ - X + n. (5) 

and t-chain, 

[a + t)^ = \2{-Tx ^t"^ + 11 + 2X)a - 2tTx + + 2r^ 
At last, after introduction of the auxiliary function (denoted by ^), 

cr r = 

we obtain the LA pair for the (4) 

^^ = -^2^2r*-A + /x, *t = 2f(r2-r^ + /x + 2A)^ + 2(/i-A)rl . (6) 



Starting out from the (6) we can find the m^KdV (which be reduced to the exponential 
CD equation by an exponential point change of variables) with it's LA pair. Continuining 
this procedure, we obtain the elliptic CD equation (m^KdV) and it is the end of this way. 
The m^KdV equation is the last equation which can be obtained by this method. 

At first, we may wonder why the limitation appears by the n=3. However, we can write 
the LA pairs for the m°KdV with n < 4 as two Riccati equations. It is incorrectly for the 
m^KdV [2]. This equation, therefore, is different from the ones found earlier and cannot 
be used by the old scheme. In the other hand, as we have mentioned in Introduction, the 
m°KdV equations (n=0,...,3), together with the Krichever-Novikov equation, exhaust all 
the integrable equations of the form ut + u^xx + fiuxx, Ux, u) — 0. So this limitation is 
something behind that. We can understand the triviality of this limitation with the help 
of the Painleve property. Let consider the m^KdV equation (4). The function r is the 
solution of both this equation and two Riccati equations (3) (inserting from the first 
equation (3) in the second one we obtain the second Riccati equation), so it is obvious 
that the function r has the Painleve property and that it is incorrectly for the solutions 
of the m^KdV with n > 3. 

In [5] the continuation of the KdV equation hierarchy in direction corresponding to 
the negative power of the spectral parameter is constructed by the following way: let the 
L-operator has the form (1), 

L = -dl + u{x,t), 

and the nonlinear equations are 

Lt^[L,A^]. 

The operators Aj^ have the form 

N 
m=— 1 

where = —1, —2, ... and are some operators. 

The first lower KdV (N=-l) equation has the form (see [5]), 

KdV_i(a) = {al + axx)t - (e'")^ = 0, 
where a is connected with u by the following way 

u = -(7xx - crl- 

Let a = iq/2, where q = q{x, t). It easy to see that 

KdV_i((T) = ^ {idx - qx) {qxt - 2 sing) - 0, (7) 

so wc have the Miura transformation between the KdV_i equation and the sine- Gordon 
equation. 

LA-pair for the (7) has the form, 

= - I + A^) ^, A = ^e^' - f ^) , (8) 

Starting out from the (8) we can find the LA-pair for the sine-Gordone equation. To do 
this we introduce new function ip 

= ^ (V'x - ^q.i^ . (9) 



therefore 

V^x = - ^q^2p, 2pt = ^e'V, = {iQxt - e'") ip. (10) 

Substituting 2 sing in place of q^t (in the eq. for the ipt) we can see that (9), (10) are the 
well known LA-pair for the sine-Gordon equation. 

Thus, we can use the method of dressing chains to both higher and lower KdV equa- 
tions. The sine-Gordon, KdV, and CD equations can be constructed starting out from 
the KdV equation via dressing chains of discrete symmetries. There are simply dif- 
ferent representations of the same equation. Prom this this standpoint, the difference 
between sine-Gorgon, KdV, mKdV and CD equations is similar to the difference between 
the Maxwell equations in different gauges! 

Surprisingly, much of the analysis of the simple system carries over directly to the 
multidimensional case. The KP equations are given by 

ut + 6uu^ + Uxxx = -3q;^ Vy, = Uy, 

where u = u{x, y,t), v = v{x, y, t), = ±1. The LA pair for the KP equations is: 

Oill^y + 'ipxx + UiIj ^0, ljjt + ^i^xxx + ^Ul/jj; + 3(Wa; - av)ljj = 0. (11) 

Setting / = log0, r = fx, where is a partial solution of (11), we can write LA pair (11) 
as 

OiTy + (tx + r'^ + u)^^0 

(12) 

n + i^Txx + 12rr^ + + 6ut + - 3av)^ = 0. 
It follows from the first equation in (12) that 

u = -Tx - - aF, Fx = Ty. (13) 

Substituting (12) into the second equation in (13) and making some simple transforma- 
tions, we obtain the well-known mKP equation [6] 

Tt - Gt'^Tx + Txxx = 2>a{2TxF - aFy), F^-Ty^ 0. (14) 

We now construct dressing chains of discrete symmetries. The KP equation admits the 
Darboux transformation [1] 

1pl=1p-T1p, Ui = U + 2Tx, Vi=V + 2Ty. (15) 

Setting s = log'?/^i, a — s^, we see that a satisfies the system of equations obtained 
from (12) by replacing u ui and v vi. Comparing these equations with (12) and 
eliminating the potentials u and we obtain 

a{s-f)y + {s + f)xx + sl-fl^Q, 

{s - f\ + [2{2s + fU + 64 - ZfX + 44 + 6 (/.. - fl - afy) s^- (16) 
6a(/., - /./,) + 2/1 = 0. 



The first equation in (16) determines the y- chain and the second determines the t-chain, 
denoted by C^^^ and Cl^^ respectively. Later (in Sec. 3) we will generalize the Borisov- 



Zykov approach for the case of the KP equation, and we will find that the results of 
Borisov-Zykov carry over with only one important changes. 

As we seen, the Painleve property is connected with the dressing chains. The aproach 
of using dressing chain (5) to construct the P4 equation was proposed in [7]. It is inter- 
esting to apply this method to the chain (16). 

To do this it is helpful to incert the periodic condition. Well known that the periodic 
condion for the (5) led to the P4 equations. Let consider the dressing chain (16). 
Setting f — fn, s — fn+i we use the condition of periodic 

fn+N ^ fn + C{y). (17) 

It easy to see that 

Un+N = tin + ac', 

SO the condition (17) get us some generahzation of the harmonic oscillator. We choose 
N — 3 and c = —2y/a. Introducing new fields Qm n = 1, 2, 3 

/i = ^ (^1 - ^2 + ^3 + c) , f^ = \ idi + 92 - 93 -0), /s = ^ (-^1 + 92 + 93 + C), 
we can now rewrite the equations for the by the following way 

a {92y - 9s y - Cy) + {qb^ - 93 9ix + 9ixx^ 0> 

Oi {-9ly + 93y + Cy) + {gs^c - 91 x) 92x + 92 XX = 0, 
a{9ly- 92y - Cy) + {9lx- 92x) 93 X + 93 xx = 0. 

Excluding g^, and using the compatibility condition Dyg2,xx — D'^92,y we get the nonlinear 
equation which can be written via some transformation as 

Zxx^]:— + + 4a;^' + 2(x2 - 2)z + ^+ 
2. Z Z Z 

(18) 

H — 6aqz H — Dy yz + 2xz — saqzj , Zy — q^- 

where z = gi^ and q = g^y. The one-dimensional limit (where Dy = g = 0) of Eq. (18) is 
the P4 equation, so the (18) is the two-dimensional generalization of the P4 equation and 
it should be called the 2D-P4 equations. 

The price we pay for the D—2 is that locahty in x is lost. That is in the usual run of 
things for the multidimensional integrable systems. 



3 Conjugate chains for the KP equations 

The correspondence between the KdV dressing chains formalism that we develop in the 
previous section and the KP dressing chains formalism is quite remarkable. We find that 
almost the entire KdV formalism can be imported into the KP chains. 



Chains (16) involve two LA pairs for Eqs. (14), which arise after introducting the 
auxihary function (denoted by ^) according to one of two rules, 



* = * = « + (19) 

Using the first rule, we obtain 

a% + (2/ + *),, + 2/.*, + = 0, 

% + 2 (2*,, + 6/,*, + + 3^)^ + Q0^^ + 4 (3/ + *)^ = 0, (20) 

= fxx + fl - afy 
from (16). The compatibility condition for these equations is 

{ft - 2fl + U^^)^ = 3a {2UJy - afyy) , (21) 

which is merely another form of the mKP equation (for r = and F = fy). 

This process can be repeated. In particular, the system (20) contains a new nonlinear 
equation obtained by eliminating the potential /. For this purpose, we linearize the first 
equation in (20) by the substituting 

/ = -^(* + «0, (22) 

where ^ = ^{x,y,t). Inserting (22) in the second equation in (20), differentiating it with 
respect to x, and introducing S = "ifx, we obtain the system 



St + 5... - IS'S, = -3a' [(le^ + QS+ {lex + 

Sy = {^XX + ixS)x , 



(23) 



Although Eqs. (23) look hke a two-dimensional generahzation of the mKdV equation, the 
one-dimensional hmit (where dy — 0) of (20) is 

(where g{x,t) = C,x = exp(— \I')) rather than the mKdV equation. After an exponen- 
tial point change of variable, this equation reduced to the exponential CD equation. 
Therefore, Eq. (23) should be called the two-dimensional CD equation rather than the 
two-dimensional mKdV equation. 

In [2], another representation for * was used in the derivation of the CD equation, 
namely, the representation determined by the second formula in (19). We have used this 
representation in [4]. 

In addition to the usual DT, the KP equations admit the so-called binary Darboux 
transformations [1], which, as we now show, allow us to construct new KP dressing chains. 
It is obvious that the KP equation admits the LA pair 

-c^Xy + Xxx + ux = 0, Xt + 4:Xxxx + QuXx + Hux + av)x = 0, (24) 



which can be obtained from (11) by simple replacing a — > —a. Let 

1 / 3 \ 

dQ{x, ^) = xHx + - iXx'ip - xi^x) dy + 4: Uj^Xx - Xxx^ - x4^xx - ^uxi^j dt 

Q{x,i')^IdQix,^)- 

It is easy to see that the one-form dQ{x, ip) is closed if ip and x solutions of Eqs. (11) 
and (24). Pair (24) also admits the Darboux transformation 

X-i = X - px, u_i^u + 2p^, v_i^v + 2py, (25) 

where p = (logx)a;, and x is a partial solution of (24). Applying DT (15), we have 

A + BQ{x,i^) 
Xi = ; , 

where A and B are arbitrary constants. Now using (25), we obtain 

= u + 2 [log(A + BQix, , ^^1,-1 = v + 2 [log(A + BQix, ■ 

We set a = Xi,x/Xi- This function satisfies the equation obtained from (12) by replacing 
u Ml, V ^ vi and a — >■ —a. In other words, this yield new integrable y- and t-chains 
(C^~) and "*), which can be called conjugate to the C^"*"^ chains discussed above: 

-a{s + f)y + {s + fU + sl- fl = 0, 

{s + /), + [4(s + fU + 64 - + 44 + 6 (/,, - fl - afy) - Gaf^fy - 2fl = 0, 

where we still have r — and a — s^- Comparing the expressions for C^~^ and C^'^^ 
from (16), wc see that these chains do not reduce to each other, thereby justifying the 
term "conjugate". This implies that the chains C^^^ lead to new integrable equations 
and also to the corresponding LA pairs. To construct these, we should again use two 
possibilities to define ^, see Eq. (19). In the present case, it is convenient to choose the 
second formula. As the result, we obtain yet another LA pair for the mKP equation (21): 

-a^y + - 2/^*^ + = 0, 

% + 2 (2*,, - S/.^r, + 3*2)^ + 6 (/2 - afy) - 6/,*,, + 4 (* - 3/), = 0. 

Linearizing the first equation by the substitution that differs from the (22) only by the 
sign of the right-hand side and defined S* = ^^a,, we get 

St + ^S^xx - IS'S, = 3a + 2^,S, + a {^y - ^C) S] , 

(26) 

aSy = {Sj: - a^xS)^ . 

Equations (26) are similar to (23) and also reduce to the CD equation in the one- 
dimensional limit (for = 1). At the same time, as follows from the method of their 
construction, Eqs. (26) do not reduce to system (23) via a 1:1 change of dependent 
variables (by a gauge transformation); therefore, Eqs. (26) can be called conjugate to the 
Eqs. (23). 



4 Dressing chains for the NLS equation 

The results in last section indicates that generahzing the dressing-chain method to con- 
siderably more complicated NLS and DS equations can be very fruitful. In this section we 
use the dressing chains of discrete symmetries to proliferate the NLS equation. Recently, 
Shabat has studied such chains for the Zakharov-Shabat spectral problem in some unusual 
gauge [8]. We have another aim: to construct higher mNLS equations. We can do it in 
two gauges, each with its own advantages and disadvantages: 

(1) Standart symmetric Zakharov-Shabat gauge. In this gauge we have very clear con- 
nections between NLS and m'^NLS equations. The disadvantage of this gauge is a very 
rapidly growing ammount of calculations. For the n > 2 we obtain very cumbersome 
equations. 

(2) New Shabat gauge [8] . The advantage of this gauge is that equations has more compact 
form. However, the connection between real NLS and real mNLS equations is not obvi- 
ous in this gauge. The reconstruction of the m°NLS for a n > 2 must be checked tediously. 

Here, we choose the standard gauge for the spectral problem and therefore the LA pair 
for the NLS equation is 

= -2ia3*A^ + 2iU^A + V^, = -iaa^A + it/*, (27) 

where 

u-ill)' "-(i-O- ^-^^i'^-^')- 

and A is arbitrary constant matrix. The compatibility condition get us the NLS equations 

iut + Uxx + 2it^i> = 0, —ivt -\- Vxx + 'iv^u — 0. 

Let $ is a partial solution of (27) with the A = diag((A + /u)/2, (A — /u)/2), * is the same 
for the Ai = diag((Ai + /ii)/2, (Ai — A''i)/2) and r = $A$^^. This matrix function is the 
solution of the system 

Tx^i[T,a^]T + i[U,T], Tt^2TxT+[V,T]. (28) 

The Darboux transformation has the form 

^^^i = ^Ai-r*, U ^Ui = U +[T,az]. (29) 

There are two discrete symmetries of the Zakharov-Shabat spectral problem: S- and T- 
symmetries. The transformation (29) is S"^ symmetry. We must use squared S-symmetry 
because an elementary S-symmetries are unsufficient to construct exact solutions of the 
NLS. 

Using (29) we obtain 

u ^ Ui — u — 2b, V ^ V\ — V -\- 2c, (30) 

where 



(31) 



Substituting (31) into (28) we get 

ttx = —dx = —i{2bc — UC + vb), — —i{2bd + (a — d)u), 

Cx — i{2ac + (a — d)v) 

and 

a* = {a^)x + 2fexC - bvx - cu^, dt = (ci^)^ + 2bcx + bv^ + cu^ 

bt = 2{ibuv + Gxb + bxd) + (a — d)ux, Ct = 2{—icuv + acx + cdx) + (a — d)vx 
Calculating the determinant and trace of matrix r we get 



(32) 



(33) 



— /x^ 

ad — be— — , a + d — X. (34) 

Using (34) we get 

a = ^ (^A ± sJJ^fi^^c^ . (35) 

Eliminating u and v from the last two Eqs. (32) and putting everything together we find 
mNLS (m^NLS) equation 

(//^ - 46c) {ibt + bxx - 26^c) + 2A (Ac + 2icx) 6^ + 2{bxC + 2bcx)bx = (ju^ - 46c) /3(6, c) = 0, 

(^2 _ 4?,c) (-ict + Cxx - 2c%) + 2A {Xb - 2ibx) (? + 2(6c^ + 2bxC)Cx = (/x^ - 46c) 7(6, c) = 0. 

(36) 

The designeshions (3{b,c) and 7(6, c) will be useful (see (41), (42)). 

As we have mentioned above, the NLS equation has a broader set of symmetries. In 
addition to the Darboux transformations, the discrete symmetries include the Schlesinger 
transformations or T-symmetries. The T-symmetries produce an explicity invertible 
Backlund auto-transformations for the NLS equation 

1 

U ^ Ui = Uxx + U^V -■! V ^ Vi — 

u u 

and 

1 v'^ 

U — > U-i — -, V ^ V-1 — Vxx + V^U 

V V 

where 

{Ui)^l = (W-l)l = u, {vi)-i = iv-i)i = V. 

These symmetries are related to the Toda chain and some of its generahzations [9] . 

The mNLS equation (36) has the similar property. We have very cumbersome Schlesinger 
transformations for the (36) so we show ones when X — /i — 0. In this case, the (36) has 
an elegant form 

ib, + bxx-2b'c-^-^-l-^f^0, 
c 2 6 

(37) 

-ic, + c..-26c^-^-i^ = 0. 

6 2 c 



The explicity invertible Backlund auto-transformations for the (37) are 

1 (2cbl - 2cbb^^ + cJJ + AcH^)^ 4cb^ 
b^bi = - -TjT^ — T , c^ci- 



(38) 



4 c^b{bl-Acb^) ' ' bl-Acb^' 

Abc^ _ 1 {2bcl - 2bcc^^ + b^c^c + Ab^c^f 

It is easy to verify that 

(6i)_i = (6_i)i = 6, (ci)_i = (c_i)i = c. 

To obtain dressing chains of discrete symmetries S'^ we construct new matrix function 
Ti = ^lAi^f^ with ^1 from the (29). Its elements 04, bi, Ci and di are solutions of 
(32)-(35) by replacing A — > Ai, /i. — > Hi, u ^ Ui and v ^ Vi (see (30)). Eliminating 
potentials u and v, we obtain our chains 

(A - 2a)6i,^ - (Ai - 2ai)b^ + 2i(Ai - ai)(A - 2a)6i - 2i(Ai - 2ai)ab = 0, 
(A — 2a)ci^x ~ {^1 ~ 2ai)ca; — 2i(A — 2a)aiCi + 2i(Ai — 2ai)(A — a)c — 0, 

bi,tb - bibt + 2 [(A - a)bib^ - (Ai - ai)6i,^6] + (Ai - 2ai) (fe^) + Xifei + M^b = 0, 

c^j^ + X2C1 + M2C = 0, 

where 

Ki = 4iC(2, 2)62 + 2[(a - ai)^ + 2iB{-2, -2)c\b - (A - 2a)B^{-2, -2), 
ir2 = -4iS(0, 2y - 2[(a - oi)^ + 2iC(0, -2)6]c - (A - 2a)C^(0, -2), 
Ml = (Ai - 2ai)S,(-2, -2), M2 = (Ai - 2ai)a(0, -2), 

^) = A^2^ ' ^) = A^2^ ' 

and fields a, d (oi, di) are expressed in term of b, c {bi, ci) via (34)-(35). 

Chains (38) involve LA pair for the mNLS (36), which arise after introducing the 
auxiliary fields 

V^=-r, 0=—, A = ai. 
b c 

After some calculations we obtain 

v. = (2^{A - AO - I) « + '-^^^^B{0, 2), 

A — 2y4 

ijj^ = 2(Ai - A)^^ - Aibc(l)tp^ - 2i{bC{2, 2)i) + cB(0, 2)0)^/' + + ^— 5^(0, 2), 

A — 2A 

(pt = 2^0^ + 4ibcip(f)'^ + 2i{bC{2, 2)iP + cB{0, 2)0)0 + Q(p - — C^(2, 2), 

c 

= ~i[2bc(pij + bC{2, 2)^ + cB{0, 2)0], 
A = -4:iXibc(j)ip - b[2iAC{2, 2) + C^(2, 2)]ip + c[2i(A - Ai)B(0, 2) + B^(0, 2)]0, 

(39) 



where 

[26^(Ai -A)- tb^^]{i^^ - 46c) + 2 (bc)^ [2b{X - a) - tb^ 



P 



Q 



b _ ibc) 

{2c^A + ic^^){tJ? - Abe) + 2 ibc)^ {2ca + ic^) 



c {ii^ - Abe) 

The mNLS equation (36) arise from the compatibihty condition of (39): 

= (V'*). , (0.), = (00. , {A^)t = (A), . (40) 

From the first two equations (40) we get two nonhnear equations 

/3i(6,c) = 0, 7i(6,c) = 0. 
The connection between these equations and the mNLS (36) is given by 

ft(6,c)=t(^) . ,.(t.c) = c(lM) . (41) 
so we have not (36) but 

p{b, c) = mb, -f{b, c) = m'c, (42) 

where m m' are arbitrary constants. Substituting bt and q (which are expressed from 
(42)) into the third equation of the system (40) we obtain m' — m. Thus, using the gauge 
transformation 

b e"'™*6, c ^ e*™*c, 
we can reduce the (42) to mNLS (36) 

(5{b, c)-mb^ e-*"^*/3(6, c) = 0, 7(6, c) - mc ^ e'^^^^-fib, c) = 0. 

This completes the proof that the system (39) is the LA pair of mNLS. We are giving two 
formulas (for convenience of reader) which are useful to check of this LA pair: 

^^""^^ at^ Xa^ + b^c-bc^-bC^iO,-2) + cB^{-2,-2), 



A -2a' 



We will shortly see the advantage of carefully working out the details of the NLS 
dressing chains. We will find that almost all of this formahsm carried over directly into 
the more complicated DS dressing chains. Formulaes for the DS dressing chains are very 
bulki but ones has the same structure. 

5 Conjugate chains of the DS equations 

The DS equations has the form 

lUt + u^x + -^Uyy -u v + qu^O, -ivt + + -^Vyy nV u + qv ^0, 

Qyy - a^Qxx = -4 (uv)^^ , 

where u = u{x,y,t), v = v{x,y,t), q = q{x,y,t). We have the DS-1 system if a = 1 and 
DS-II, ii a — i. Under the reduction v — ±u, one obtains a known model that describes 



the propogation of a small amplitude wave packet that is quasi-one-dimensional and quasi- 
monochromatic over the surface of a nonviscous curl-free liquid [10]. The quantity u is 
the envelope of the velocity potential, while q describes the nonlocal flow generated by the 
wave packet. A different application of (43) is related to the dynamics of plasma waves 
[11]. We woun't use this reduction restriction in this section. 
LA "pair" for the (43) has the form 



i/jy = ai/j^ + u(f), (f)y = -a(f)x + viIj, 



2i 

ipt = 2iip^^ H u(l)^ + ( ^ 

a 

2i 

(t)t = -2z0^3; H vil)^ + - 

a 

















a 





—F — F 



% \ % 

—uv + — {au^ + Uy) (f), 



+ —uv I + ^ {ciVx - Vy) ip, 



(44) 



a 



where q = —iF^. 

Let {ipi, 4>i; ip2, 4>2] ipj 4>} solutions of (44) for the same u, v and F. DT is given 

by 



ip ^ ipi = - a^p - b(j), (p ^ 4>i = 4>x - cip - d4>, 

u ^ Ui — u + 2ab, V ^ Vi — V — 2ac, F ^ Fi — F + Ai{a + d), 



where 



i>l,x<l>2 - i>2,x<l>l , i>2,xi>l - V'l,xV'2 01,x02 - 02,101 

— ■ 1 ■ , 0= — ■ 1 ■ , c— — ■ : ' . 



d 



(p2,x'4'l - 4>l,x'4'2 



A = -0102 -'0201- 



The quantities a, b, c d are solutions of the system 

Gy = a{ax + 26c) + uc — vb, by = a{bx + 2bd) + {d — a)u + Ux, 

Cy — —a{cx + 2ac) -\- {a — d)v -\- v^, dy — —a{dx + 2bc) — uc + vb, 



(45) 



(46) 



at — 2i 

i 



ax + 2bc + a'^ + ^] + 2{a + d)bc 



4 



+ 7? [K + cuy- {uv)J[ + 



+ - 



a 



cUx — bvx + 2{cdu — abv + {uc)x) 



bt = 2% 
i 

+- 



{bx + 2bd)^ + 2 ( + 6c + + I ) 6 



H — J [uxy + {d- a)uy - 2buv] + 



a 



a 



Uxx + {^d — a)ux + 2(dx + d — ad + bc)u — 2b v 



Ct = -2i 

i 

+ - 
a 



dx + be + + -j c - — [vxy + (a - d)vy - 2cuv] + 
Vxx + (3a — d)vx + 2{ax + — ad + bc)v — 2c^u , 



(47) 



dt = -2i 



dx + 2bc + d'^ + '^] + 2{a + d)bc 



— [bVy + CUy- {UV)^] 



a 



I 

a 



q 

cUx — bVx + 2(cdu — abv — {bv)x) — ^ 



These equations can be obtained from the (44). Eqs. (47) are similar to the (13)-(14). 
The LA pair (44) and Eqs. (47) are simply two different reprersentations of the same 
equations. To construct mDS equations we must, first, express fields a and d via one 
function 9 — 9{x, y, t) defined by the relations 

2q; ' 2a ^ ^ 

It is easy to check the truth of this representation by the substitution (48) into the first 
and fourth equations (47). Second, we must exclude potentials u and v from the first four 
equations (47). Introdusing new fields ^ — ^{x,y,t), X — 9^, Y — 9y 

q = -2X, _ - 46c + ^ (e - 

and putting everything together we find mDS (mDS-I if a = 1 and mDS-11 if a = i) 
equations as the system of five equations for the five functions b. c. ^, X and Y: 

ibt + 2{bX + b^)^ + i [{UX + U^)^ + 2b {Uc - Vb) - 2Y {bX + b^)] + 

+ A [Uxy - U{Xy + XY) - 2{U,Y + bUV) + ib\ + \y {UY - Uy) = 0, 
a^ a"^ 

-ict + 2 (cX + c,), - ^ [{VX + V^)^ - 2c {Uc - Vb) - 2Y {cX + c,)] + 

+\ [Vxy - v{Xy + XY) - 2(y,Y + cUV) + Cc] - \y (VY - Vy) = 0, 



iXt + 4 [(bc)^ + 2Xbc] - - [X(Vb - Uc) + Achy - 2(cU)^ + bV^ - cU^ - YX^ 

2 ^ 1 

-- [{UV), - Vby - Ucy + YibV + cU)] - (y'-c)^0, 

iYt + A{byC - bCy) + ^ [2{UyC " U Cy + Vyb " Vby) + Q = 0, 

~Yx — 0, 

where 

U= „ ^\ V 



(49) 



2a{2bcY + a{bc^ - 6^c)) ' 2a{2bcY + a{bc^ - b^c)) ' 

Gi = a^%X^ - bX^^ + 2b{2bcX + b^c - bc^)] + a^[2{bcy - byc) - {X^ + 46c)r]6+ 

+ a{bXyy - b^Yy) + bYYy, 

G2 = a^lc^Xa; - cX^a: + 2c{2bcX - b^c + 6c^)] + a^[2{bcy - byc) + {X^ + Abc)Y]c+ 

+ a{cXyy - CxYy) " CYYy. 

As we have mentioned above, the NLS equation has two types of discrete symmetries: 
Darboux and Schlesinger transformations. In just the same way, the discrete symmetries 
of the DS equations include T-symmetry. This transformations are given by ([12], [13]) 

u ^ ui^ u{uv ^ a^ilogu)^:^- {\ogu)y^ , v^Vi^-, q ^ qi ^ A{\ogu)^^ 
u^u^i^-, V ^v^i^v (uv + a'^ {log v)^^- {log v)yy^ , q ^ q^i ^ q + A{logv), 



The similar transformations for the (49) are given by 



h h - ^ 

' aU{U + 2aby ^ U{U + 2aby 

V ^ - ^ , o UK-Uxb ^ _ -U^Y + 2a[UU, - hUy + aU% + hX)] 
X^X,-X + 2ajj^^-^^, Y^Y,- [/([/ + 2ah) 

, , a'U{b,. + {bX%) + a[b{UY^ - U,y) + U{U,, - b,Y)] - UU,Y 
e-6-e + 4a + 2ab) ' 

and 

C — >• C_i = r, — >• 0_i 



ay (V^ - 2q;c) ' y(\/-2Q;c)' 



V^(\^ - 2ac) ' y (y - 2ac) 

a^V{c,, + (cXy) - a[c{VY, - K,) + V{V,, - c,Y)] - VV,Y 



V{V - 2ac) 
where 

Ml = [U%{2X, - X^) + 2U{{U,b - Ub,)X + {U^)^ b, + A{Ubfc + (t/J - 2UU,,) b] + 
+q;2C/ [C/(C/X^ + 26y^) + (2C/j^6 - UU^)X - 2Ub^Y + 4C/2c6 - 2U^yb - UU^^ + 2b^Uy + 2C/2] + 
+a [U\Y, + + C/c) + C/2(6y2 _ 2U,Y - U,y) + Uy ((C/^), - 6C/^)] + U''Y{UY - Uy), 

M_i = a' [V'ciX' - 2X,) - 2V{y,c - Vc,)X - {V^)^ - Aiyc)% - {V^ - 2VV^,) c] + 
+aV [V{VX:, + 2cFa:) + (2l^j,c - VV^)X - 2Vca;F + AV%c - 2V^yC - VV^^ + 2c:,Vy + 21/^] _ 
-a [V\Y, + XY + Vb) + VicY-" - 2V,Y - V^y) + Vy {{V^), - cVy)] + VY{VY - Vy). 

It is possible to verify that X±i^y = Y±i^x, (&i)-i = (&-i)i = b and so on. 

To obtain dressing chains of discrete symmetries we must exclude potentials u, v, 
q from the two systems: (47) and and the system of equations obtained from (47) by 
replacing a — > ai, 6 ^ 6i, c — > Ci, d — > di, g — > = 5 + 4(a + d)^, u ^ Ui, v ^ Vi (it is 
necessary to use (45)). As a result, we have y-chains 

[2bc{d -a) + bxC - bc^] {ai^y - aai^^) + {bci - bic){ay - aax)x + {bci - bic){bcy - byc) + 
+a[{bic + 2biCi + bci){bcx - b^c) + 4bc{aci{b + 61) - 6id(c + Ci))] + 
+[(6ci + bic){a - d) + bic^ - bxCi\{ay - aa^) = 0, 

(6ci - bic) [{bi - b)y - a{b + bi)^] + (a - 04 - d + di) [a{ai^xb - aj)i) + aybi - 01,^6] + 
+2a[{b + bi){abci + cbidi) - bbi{c + ci)(ai + d)] = 0, 

{bic - 6ci)[(ci - c)y + a(c + ci)^] + (a - ai - d + di)[a{ai^xC - a^ci) + QyCi - ai^yc]- 
—2a[{c + Ci){aiCib + cdbi) — cci{b + bi){a + di)] — 0, 

(ai + di)y - a{ai - di)^ = 

(50) 



and t-chains: 

Q;^(ai — a)t + i[2a^((a^ — — 2hiC\ — ai^^ — dx)x + 2((a + d)hc — (6 + h\)cid-i — (c + Ci)ai6i — 

-{hci),j) - b^ci - biCj.) + a((c - Sci)?/^ + (& + ^OK- + 2(6iCj^ - byCi - (a + cZ)^j,+ 
+ aibi - ab)V - (ci,^ + cidi - cci)C/)) + (6 - bi)Vy + (c - Ci)C/j/] = 0, 

a\bbi,t - btbi) + i[2a2(262(ai(ii - bid) - 2bbl{ci + c) + b^^xbi - 6i,xx& - 26o?i,^(6 + 
+2(6^d6i - + b[aibx - 2ai,^6i - b^^ - Sb^di + 2{bi{d^ - be) - dl{b + 6i))])+ 

+Q;((fei - &)f/x.x- + (3(&irf - + ai6 - abi)Ux + 2[6i&\/(6 + bi) + 66j/(ai - di) - bb^y^ 
+U{bidx — bdi^x + bid^ — bd\ — bbi{c-\- Ci) + ba^di — b^ad)]) + Uxy{bi — b)-\- 
+Uy{b{ai - di) - bi{a - d))] = 0, 

oP'icci^t - CtCi) + i[2a^{2c^{bici - aidi) + 2cc\{bi + 6) + ci^^^c - cic^^ + 2cai^x{c + Ci)+ 
+2(caiCi,a; - ciac^) + c[2cidi^x - dic^ + c^x + ^c^ai + 2(ci(6c - a^) + ai(c + ci))])+ 
+a((ci - c)\4j, + (3(aci - aic) + cdi - Ci(i)\4 + 2[cic[/(c + Ci) + cCy{di - ai) - cCxy+ 
+V{ciax — cai^x + cio^ — ca\ — cci(6 + bi) + caidi — ciad)]) + Vxy{c — ci) + 
+Vy{c{ai - di) - ci(a - d))] = 0, 

— d)t + ?'[2a^(((i^ — ci^ + 26iCi + + (ii^x;)^ — 2((a + d)bc — (6 + 6i)ci(ii — (c + Ci)ai6i — 
-(&ic)x) + &xCi + biCx) + - 36i)K + (c + ci)?7a, + 2(ci&y - biCy - (a + 
+(c^ + cidi - c(i)t/ - + ai6i - a6)l^)) + (ci - c)t/j/ + (6i - b)Vy\ = 0. 

(51) 

These formulas have awful sight! Unfortunately, I don't know how to write them in more 
compact form in our gauge. (51) has the best form which I could imagine. 

Although this equations are long and difficult, the end result is quite simple. We can 
see that all of (1 + 1) NLS formalism carried over directly into the (2 + 1) DS dressing 
chains. In particular, the chains (51) involve LA pair for the mDS (49). We do not give 
it here. 

As noted in Sec. 3, the KP equations admit two types of chains, which we call 
conjugate. It is a new result, which is characteristic precisely of multidimensional systems. 
Of course, it is truth for the DS equations. 

To construct the chains which are conjugate to (50)-(51), we must introduce the new 
LA pair for the DS equations {p = p{x, y,t), f = f{x, y,t)): 

Py = aPx - Vf: fy = -a fx " Up, 

2i /Ifl \ i 

Pt = -2tpxx + -vfx - i^- ^-Fy + F^J - —uv )P + -^ {avx + Vy) /, (52) 

2i /Irl ^ i \ i 

ft = 2ifxx + —upx - - -Fy - Fx + -^uv / + ^ {aux - Uy) p, 



two matrix functions 



and one-form 

dVt = ^"i/dx + a^aa'i/dy + 2i ^$0-3 ^^o, - ^xC^a"^ + -$(73^^ , fl = J dfl. 



This one -form is closed if V'1,2, <t>i,2 , Pi, 2 and /i^2 are solutions of Eqs. (44) and (52) for 
the same u, v F. It is easy to see that the quantities 

. Pl,xf2 - P2,xh n h,xf2 " f2,xfl ^ P2,xPl " Pl,xP2 

^ — 1 5 ^ — ^ 1 ^ — ^ 5 

A A A 

/2,a;Pl — h,xP2 \ r r 

D^— ^-^ — , A = pi/2-p2/i- 

are solutions of the y-system 

Ay = a{A^ + 2EC) + -uC - vB, By = -a{B^ + 2AB) + {D - A)u - 

(53) 

Cy - a{C, + 2CD) + {A- D)v - v,, Dy = -Ay + a{A - D),. 

Here, we don't need t-systcm (four equations for the At, Bt, Ct and Dt). 

The Darboux transformation (45) (for the system (44)) give us the transformation 
rule for the functions A, B, C D 

A ^ ai ^ All - a, -B ^ 61 = A12 - 6, C ^ ci = A21 - c, D ^ di^ A22- d, 

where a, 6, c and d are defined in (46), and Aj^ are elements of the matrix A — 
i,k = 1,2. Transforming the system (53) and exluding potentials u and v from this system 
and from the first four equations (47) we get y-chains which are conjugate to (50). There 
are four equations, two of them arc equaivalent to the first and fourth equations (50). 
The rest two equations are given by 

{bci - bic)[{b + bi)y + a{b + bi)^] + {a + ai-d- di)[a{a^bi - Ui^^b) + ai^yb - Uybi]- 
—2a[{b + bi) {abci + aibic) — bbi{c + Ci){d + di)] = 0, 

{bic - bci) [(c + ci)y - a{c + Ci) J + {a + ai-d-di) [a{a^Ci - ai^^c) + ai^yC - ayCi\+ 

+2a[{c + ci){bcidi + bicd) — cci{b + 6i)(a + oi)] = 0. 
The conjugate t-chains can be obtained by the same way. We do not give them here. 

6 On the localized solutions of the DS equations 

Discrete symmetries are a good way to obtain exact solutions of the nonlinear integrable 
equations. In [14] we applied the DT (S-symmetry) to construct exact solutions to the 
DS-I and DS-II equations (the same to the BLP equations see in [15]). In particular, we 
have obtained the dromion solutions of the DS-I equations. In this section we use the 
T-symmetry to construct nonsingular solutions to the DS-I that fall off according to the 
exponential and/or rational low along all directions in the plan. In addition to the known 
dromion, we give several new solutions with above properties. In rest path of this section 
we present (via S-symmetrie) novel exact solution of the DS-II equations describing the 
soliton on the plane wave background. 
1. DS-I equations. 



To study the DS-I equations it is convenient to use the following change of variables: 

11m 

We are interested of locahzed solutions to the DS-I equations which move without 
shape distortions and we look for these solutions in the form 



u{t, X, y) = v{t, x,y) = U (C, ?]) exp(i^), q{t, x, y) = Q{^, ?]), 

where 

^ = X — 2at, 1] = X — Aht, 9 = ax + by — {a^ + b'^)t, 

and U{^,ri) is the real function. It was shown in [16] that a solutions to the DS-I equations 
can be found in this setting from the nonlinear Liouville equation 

9,a,$ = ^exp($), (54) 

for the function $ = log(C/^). This result was obtained via the T-symmetries (see the 
previous section). 

It is convenient to introduce the boundary conditions for the U {x, y) as 

U{x,Q)^A{x), U{0,y)^B{y), A(0) = S(0) = C, (55) 

where A{x) and B{y) are the given functions. Now we can use the solution to the Gousat 
problem for (54) given in [17], from which we find 

U(x v) = ^CA{x)Biy) 
Using T-symmetries, one can also find the explisit expression for the Q{x,y) field 

Qi^^y) = l(P'ij-X + G'i^)") - ^([(logp')f + [(logG')f )+ 

, 2 , {PG'f + {P'Gf + ^P'G\ 



where 



Pix) = ^£dpA\p), G{y) = ^ f^dqB\q). 



Expression (56) also allows us to construct localized solutions. For example, choosing 
the boundary functions in the form of two solutions of the NLS equation, 

we obtain the well-known " one-dromion" solution [18]. Other localized objects arc found 
by choosing the boundary conditions appropriately. Thus, assuming A{x) and B{y) to be 
solutions of the KdV equation, we obtain a new solution that decays exponentially in all 
directions, 

U(x y) = ^ 

(cosh(aa;) cosh{(3y)y[Aa/3 - ananh(aa;) tanh(/3y)(l - tanh^("^) )(i _ tanh^C/?^) )] ' 

(58) 



where a, a and /3 are real constants satisfying the condition d} < 9a(3, which guarantees 
the nonsingular behaviour of (58). Let us note that if if the functions A{x) and B{y) have 
maxima m and /, respectively, the solution U {x, y) determined by Eq. (56) corresponds 
to ml localized bound formations. It possible to show that these equations are different 
from the (L,M)-dromions [19] and A'"^-soliton solutions built in [20] and [21]. 

The rational localized solutions ("lumps") are obtained in the same way. let us con- 
sider, for example, the following boundary conditions: 

Ai{x) = , A2{x) 



x"^ + 0? x'^ — iaxY + a'^ /mx 

Substituting Ai and Bi into (56) and using (55), one obtains the one-lump solution, 

if = -16C anY(C" 6' 1/ X - 16 ( + ) (|/2 + 62 ) 

+ ahC'^ {SQ{x^ + a^){y'^ + b^)+yhS{x^ + a^)+xaQ{y'^ + h'^))), 

S — arctan ^— ^ , Q — arctan . 

Substituting Ai and B2 (or A2 and Si), we get the two- lamp solution, and substituting 
A2 and B21 the four- lamp solution. One can always choose the paramerers involved 
in (59) such that the solutions falling off in all directions according to rational law are 
nonsingular. Finally, choosen A{x) from (57) and Bi{y) from (59), we obtain the localized 
solution falling off according to the exponential law as a function of x and according to 
the ratianal law function of x 
2. DS-II equations. 

We couldn't find a localized nonstationary solution of the DS-II equation. We present 
only soliton on the plane wave background. We choose the initial solution of the DS-II 
equations as 

w = ^ exp {iS) , S = -{2A^ + - h^)t + ax + by, g = 0. 
where A, a and b are real constants. The solutions of the LA pair (44) 



^: = /exp I if + m) , V. = ':<2n^^ (f - 2i^] oxp (if + M 



where 



X + - ((6 - 2m)y + 2 (2{bm - - 2m^) - ap) i 



1 



+ 0L2, 



2 



M = - ((p - a)y [p(6 2m) - 4am]t) , 4/1^ (6 - 



2 



0. 



Using (45) (with m — m., a.\ — a\ and 0:2 = 0) we get nonsingular solution u\ which is 
the soliton on the plane wave background: 



\2 \2 X. 2,2 

Ux \ ^ A at >oo 



This solution is two-dimensional generalization of "exulton" solution of the NLS built in 
[!]• 



7 Conclision 



The past years of intense theoretical research have made it increasingly clear that the 
secret to integr ability most likely lies in the power of discrete symmetries. Let us sum- 
marize some of the promising features of the discrete symmetries: 

(1) It is the an extremely power method to construct exact solutions of integrable equa- 
tions. 

(2) This approach includes all "sohton miracles": finit-gap solutions, the Painleve prop- 
erty cts. 

(3) Discrete symmetries allow us to proliferate integrable equations. For example, the 
whole MKdV theory can be consider as the DT theory and the Toda-Volterra theory can 
be consider as the theory of S- and T-symmetries of the NLS equations [9] . 

(4) Discrete symmetries led to connection between integrable system and supersymmetry. 

(5) Discrete symmetries allow us to construct discrete integrable systems. 

Ideally, we would want a unified theory to unite and understand all "sohton miracles" . 
There are two ways to find this unified theory. The first path is connected with the Hirota 
bilinear difference equation [22]. This famous equation is known to provide a canonical 
integrable discretization for most important types of soliton equation. 

The second path is the theory of dressing chains. This approach allow us to proliferate 
integrable equations and, at the same time, to establish a links between known integrable 
equations (it is clear that the Miura map can be obtained from the dressing chains). 
Ideally, we want to show that all integrable equations are nothing but different forms of 
a single equation! 

Even if this direction is right, it is still a long way off. Now the theory itself often 
seems like a confused jumble of random (but usefuU) rules and random (but remarkable) 
observations. It remains to be seen how useful the dressing chains will become in the 
future. 
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